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Введение. Повышение долговечности лопаток газотурбинных двигателей (ГТД) достигается за счет применения 
специальных защитных покрытий их поверхности. Для разработки таких покрытий основной исходной 
информацией является геометрический профиль сечения лопатки. Для передачи заданного профиля сечения 
лопатки в соответствующую САО/САМ-систему или пакет инженерного анализа применяются методы 
параметрического моделирования, позволяющие автоматизировать данную операцию. Однако известные 
подходы к созданию параметрической модели профиля лопатки не лишены ряда недостатков, и общепринятой 
методики ее создания в настоящее время не существует. Целью данной работы является создание методики 
аппроксимации профиля лопаток ГТД, удобной для использования при последующем анализе условий работы 
специальных покрытий поверхности лопаток. 

Материалы и методы. При построении параметрических моделей профиля лопаток газотурбинных двигателей 
авторами использован метод, основанный на применении ортогональных полиномов Лежандра. Это позволило 
обеспечить высокую точность аппроксимации и построение непрерывного отображения для параметров 
аппроксимации профиля лопатки. Для автоматизированной обработки исходных профилей создано приложение 
на языке Рушфоп, позволяющее вычислять коэффициенты аппроксимирующих полиномов для линий контура 
лопатки, визуализировать результаты расчета и создавать по точкам аппроксимирующих функций ахЁЕ файл для 
передачи его в САО-систему. Далее геометрические модели лопаток использовали для решения задачи обтекания 
лопатки стационарным аэродинамическим потоком. Результаты решения этой задачи использованы при 
исследовании влияния профиля лопатки на ее остывание в аэродинамическом потоке. 

Результаты исследования. В качестве примера рассмотрены три варианта профилей лопаток, относящихся к 
разным типам ГТД. Показано, что для всех трех изученных профилей предложенная методика позволяет 
получать параметрические модели, обеспечивающие высокую точность построения аппроксимирующих линий, 
что подтверждается близкими к единице значениями коэффициентов детерминации. Для иллюстрации 
возможности использования полученных моделей приведены примеры решения задачи газовой динамики при 
потенциальном обтекании лопатки в стационарном аэродинамическом потоке. Методом конечных элементов 
рассчитаны распределения давлений и температур на поверхности лопатки. 

Обсуждение и заключение. Результаты вычислений показали, что предлагаемая методика аппроксимации профиля 
лопатки ГТД, основанная на использовании ортогональных многочленов, является удобным инструментом, 
позволяющим автоматизировать создание геометрической модели лопатки и проводить сравнение различных типов 
и профилей лопаток, решая соответствующие задачи газовой динамики. При этом для заданного профиля лопаток 
и условий работы ГТД можно получить распределение температур и усилий, действующих на поверхности лопатки, 
что необходимо для прогнозирования долговечности специальных покрытий. 
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Введение. Лопатки турбины являются наиболее высоконагруженными деталями газотурбинного двигателя, 
долговечность которых главным образом определяет межремонтный ресурс двигателя [1]. Одним из путей 
повышения долговечности лопаток турбин является нанесение специальных защитных покрытий на их 
поверхности [2]. Покрытие представляет собой сложную композитную конструкцию, для выбора материалов 
которой и определения оптимальных геометрических характеристик необходимо знать условия эксплуатации 
детали: распределение температур, давлений и сдвиговых напряжений на поверхности лопатки. 


Информатика, вычислительная техника и управление 


Б&р://уезшК-Чопзва.га 


'А4уапсеа Епятеетия Везеагсй (Козоу-оп-Оопн). 2024;24(1):78-87. е155М 2687—1653 


К форме лопаток, а также к выбору материалов, технологии их изготовления и специальных покрытий 
предъявляются самые высокие требования [3]. При проектировании специальных покрытий основной исходной 
информацией является геометрический профиль поверхности детали, поскольку на его основе задаются 
параметры технологического процесса нанесения покрытия. 

Среди характерных элементов профиля можно выделить: хорду — расстояние между наиболее удаленными 
точками при горизонтальном расположении лопатки, подобно крылу самолета, спинку (засНоп $14е) — верхнюю 
часть профиля, корытце (ргеззиге 514е) — нижнюю часть профиля, входные и выходные кромки [4]. Выбор формы 
профиля сечения лопаток производится как на основе экспериментальных исследований опытных установок и 
реальных двигателей [5], так и по результатам численных экспериментов на моделях [6]. В процессе численного 
моделирования возникает задача создания геометрической модели лопатки на основании заданного профиля 
сечения и передачи ее в соответствующую САБ/САМ-систему [7] или пакет инженерного анализа [8], в котором 
задаются параметры физических свойств материалов деталей, определяются граничные условия задачи и 
осуществляются вычисления. 

Традиционно при создании чертежа лопатки ее профиль описывается набором дуг окружностей, который 
затем по точкам передается в модуль создания геометрической модели и генерации сетки конечных элементов. 
Недостаток данного подхода состоит в том, что при создании каждой новой модели приходится вручную 
создавать в графическом редакторе ее профиль. В этой связи были предложены методы параметрического 
моделирования профиля сечения лопатки, позволяющие автоматизировать данную операцию. Одной из наиболее 
известных таких моделей является девятипараметрическая модель КАТО (Карл Ах!а|1 Тише Пез1епе), 
предложенная в работе [9]. Вместе с тем практика использования КАТО выявила и определенные недостатки, 
свойственные данной модели, в частности неудобство использования при оптимизации геометрии профиля, а 
также недостаточную точность. На этом основании авторами [10] была предложена модифицированная версия 
данной модели, включающая в себя болышее число параметров. В отечественной практике также активно 
используются методы параметрического моделирования профилей сечения лопаток. Так, в работах [11, 12] были 
предложены алгоритмы автоматизированного построения профиля лопаток с использованием аппроксимации 
кривыми из набора парабол и кривыми Ве7ег второго порядка. В работе [13] на основе данных алгоритмов 
предложена методика проектирования решетки профилей ГТД. Использование параметрических моделей 
профилей лопаток при решении задач оптимизации рассматривалось авторами в [14, 15]. 

Вместе с тем следует отметить, что существующие методики параметризации профиля лопаток не лишены 
ряда недостатков. Основной недостаток традиционных схем параметризации профиля лопаток, основанных на 
использовании кривых второго порядка, состоит в том, что для обеспечения точности описания приходится 
использовать различные функции для разных частей детали. При этом данные функции не образуют 
ортогональную систему, а значения параметров, определяемых аппроксимацией существующего профиля 
методом наименьших квадратов, оказываются скоррелированными. Это приводит к тому, что не удается 
построить непрерывного отображения для множества параметров при аппроксимации системой сечений 
трехмерной поверхности лопатки двойной кривизны. 

Также необходимо отметить, что универсальной методики параметрического моделирования профилей 
лопаток ГТД в настоящий момент не существует, а выбор конкретной методики зависит как от целей 
моделирования, так и от особенностей конкретного типа проектируемого двигателя. Вместе с тем формы входной 
информации могут существенно различаться. Это может быть и геометрическая модель в виде файла в формате 
одной из САП-систем, и просто набор точек из базы данных профилей типа [16]. 

Таким образом, дополнительным требованием к методике параметрического моделирования профиля лопаток, 
обусловленным назначением модели, является универсальность по отношению к формату исходных данных. 

Целью настоящей работы явилось создание методики аппроксимации профиля лопаток ГТД, свободной от 
вышеперечисленных недостатков и удобной для использования при последующем анализе условий работы 
специальных покрытий поверхности лопаток, что продемонстрировано на примере решения задачи газовой динамики 
для трех различных моделей лопаток, профиль которых параметризован с применением предлагаемой методики. 

Материалы и методы. Авторами использован вариант параметрических моделей профиля лопаток, 
основанный на использовании ортогональных полиномов Лежандра [17]. Для автоматизированной обработки 
исходных профилей было создано приложение на языке Руфоп, позволяющее вычислять коэффициенты 
аппроксимирующих полиномов для линий контура лопатки, визуализировать результаты расчета в форме 
графиков исходных точек и аппроксимационных кривых для спинки и корытца лопатки, сохранять исходные 
точки профиля, массив коэффициентов аппроксимирующих полиномов и коэффициент детерминации модели в 
базе данных и создавать по точкам аппроксимирующих функций 4хЁфайл для передачи его в САО-систему и 
пакет инженерного анализа. 
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Построенные двумерные геометрические модели лопаток использовали для решения задачи обтекания 
лопатки стационарным аэродинамическим потоком в приближении потенциального обтекания. Вычисленные в 
результате решения этой задачи распределения скоростей в потоке использовали при исследовании влияния 
профиля лопатки на ее остывание в аэродинамическом потоке. 

Методика расчета. Системы ортогональных многочленов, используемые в качестве основы предлагаемого метода, 
обеспечивают линейную независимость коэффициентов аппроксимации и лишены недостатка, приводящего 
к затруднениям при построении моделей, непрерывных по параметрам. Степень полиномов выбрана достаточно 
высокой — девятого порядка, так что набор из двух полиномов для двух частей профиля (спинки и корытца) 
составляет двадцать параметров, что является достаточным для аппроксимации весьма сложных профилей. Ниже 
описаны вид аппроксимирующей функции и методика вычисления оценок ее коэффициентов. 

Обозначим аппроксимирующую линию профиля лопатки посредством функции у(х). Выражение данной 
линии в случае ряда полиномов Лежандра девятой степени представим в виде: 


9 
У(х)= У а (х), (1) 
где а; — коэффициенты, являющиеся параметрами модели; /(х) — полиномы Лежандра, вычисляемые по формулам: 
Го = 1; 
У =х, 
То = (35? -1) 2; 
18 =(5%° 3х) 2; 
14 = (35х*—30х? +3) 8; 
15 = (63х° 70%? +15х)/8; (2) 
16 =231/16х° —315/16х“ +105/16х? — 5/16; 
11 =429/16х7 — 693/16х5 +315/16х3 —35/16х; 
1в =6435/128х8 —3 003/32хб +3 465/64х* —315/32х2+35/128; 


_12155 _ 6435 9009 1155, 315 
128х° 32х’ 64х’ 32х3 128х 


Пусть координаты исходного профиля лопатки заданы набором точек у? = у° (х;), 1=1,..., М. Подставляя 


значения х; в функции (2), получим матрицу Х размером №х 10. Коэффициенты аппроксимирующего полинома 
для заданных профилей лопаток находили методом наименьших квадратов из условия минимума суммы 
квадратов отклонений заданных значений уг и рассчитанных по уравнению регрессии (1): 


а-У т (уе -У(аь, ж)} —> тт. (3) 
Оценки коэффициентов а; вычисляли с помощью известной формулы регрессионного анализа: 
аи" (4) 
где а — обозначение для вектора оценок коэффициентов; у — для вектора точек заданного профиля уг. 


Точность аппроксимации оценивалась по коэффициенту детерминированности модели (коэффициент 
детерминации А?), который при хорошей аппроксимации должен быть близким к единице: 


в.) 


. (5) 


где У’ — среднее значение ух. 

В конструкторской практике принято входную и выходную кромки контура лопаток описывать радиусами 
окружностей. В связи с этим в данной работе рассматривались координаты точек аппроксимируемых профилей 
между дугами окружностей входной и выходной кромок. Для целей аппроксимации расположение контура 
рассматриваемой лопатки было принято горизонтальным, так, чтобы линия хорды лопатки совпадала с осью 
абсцисс декартовой системы координат. Для сохранения условия ортогональности моделей обязательным 
является нормировка размеров. В связи с этим все размеры нормировались на длину хорды лопатки, так что 
абсциссы точек профиля лежали в интервале 0,1. Это обеспечивало ортогональность функций (2) с 
коэффициентами, вычисленными по формулам (4). 
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Результаты исследования. В таблице] приведены коэффициенты аппроксимирующих полиномов для 
профилей лопаток, относящихся к разным типам ГТД: профиль лопатки компрессора С8626 [13], профиль МАЗА 
лопатки высокоэффективного ГТД высокого давления [18] и профиль лопатки-имитатора, использованной в 
работе [19] для исследования термической усталости теплозащитного покрытия, состоящего из двух 
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металлических и керамического слоев. 


Таблица 1 
Коэффициенты аппроксимирующих полиномов для профилей лопаток 
из разных источников 
Коэффициенты Ве р о 

спинка корытце спинка корытце спинка корытце 
ао 1424250 —123,0993 —1275,6471 717,8336 —1043,6707 527,0915 
а —374,3446 324,9875 3328,6126 —1868,9670 2719,5055 —1369,9895 
92 477,9332 —417,7186 —4181,1239 2338,1048 —3406,1466 1707,9198 
аз —444,4554 392,7204 3794,0526 —2108,4202 3077,0694 —1531,8025 
@4 324,6819 —292,0489 —2684,7257 1479,0172 —2164,2992 1066,1471 
а5 —189,2284 174,3673 1502,4170 —818,4219 1201,7203 —583,7501 
46 86,8260 —82,4103 —654,4877 351,5869 —518,2516 2473087 
а7 —30,0775 29,5457 211,8930 —111,9355 165,6344 —7 7,2958 
аз 7,1533 —7,3024 —46,0284 23,8450 —35,3931 16,0493 
99 —0,9024 0,9652 5,1021 —2,5861 3,8411 —1,6752 
В? 0,9999 0,9993 0,9896 0,9990 0,9982 0,9994 


На рис. 1 приведены исходные точки данных профилей и соответствующие аппроксимирующие линии 


0,0 


„я 


0,0 


в) 


тети 
и м нь... ” 
0,4 0,6 0,8 х 
а) 
— = 
2% > & 
в 
-#— 
‚вай Е ть. . Ч 
0,4 0,6 0,8 х 
6) 
в. 
\. -. 
в. 
——_“^ 
5 5, | 
0,4 0,6 0,8 х 


Рис. 1. Исходные точки и аппроксимирующие линии профилей лопаток, представленных в таблице 1, 


источники профилей лопаток: а 


[13]; 


6— [18]; в — [19] 
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Как было указано выше, целью параметризации в настоящей работе является аппроксимация профиля 
непрерывной функцией с последующим моделированием распределения внешних параметров по поверхности 
лопатки. Для иллюстрации возможности решения данной задачи рассмотрим простую двумерную модель 
лопатки в стационарном аэродинамическом потоке, в приближении потенциального обтекания. В этом случае 
распределение скоростей потока у(х, у) описывается уравнением: 


я 
Ох 
у= д | (6) 
ду 
где функция \у находится из решения уравнения Лапласа: 
Д\у = 0, (7) 


которое в данном случае отвечает условию несжимаемости (У -у= 0) и безвихревому потоку (У ху=0). 

Граничным условием для функции \у будет являться равенство нулю скорости в нормальном направлении 
поверхности лопатки, что означает постоянство \у на этой поверхности. 

Обозначим линию профиля поверхности лопатки посредством 5 и будем искать решение уравнения (7) в 
области ©, ограниченной внешней частью по отношению к 5 и внутренней частью по отношению к окружности С 
достаточно большого диаметра, по сравнению с длиной хорды лопатки, так что на поверхности С поток является 
однородным. То есть на границе области 009 = 5/С заданы условия: 


—_ — ЕВ 
[$ =0, ч[с = У» Х", (8) 
где у» — вектор скорости однородного потока. 


После вычисления скоростей можно вычислить распределение давлений по уравнению Бернулли: 


У 

ру 
=В-Р 9 
р - (9) 


2 
2 У 
где р — плотность газа; В — постоянная Бернулли, которую можно положить равной р.+ =>”, где 
Ро — давление на бесконечном удалении от лопатки. 
Вариационная постановка задачи (7) выглядит следующим образом: 


[Учун=0, УшеЕ Но (0). (10) 
о 
Решение данной задачи осуществлялось методом конечных элементов. Для сравнения расчет был проведен 
для трех вариантов профилей лопаток, приведенных в таблице 1, с одной и той же длиной хорды /. 
Полученные результаты использовались далее при решении задачи остывания лопатки в аэродинамическом 
потоке. Такая задача имеет смысл при термоциклических испытаниях покрытий на лопатках-имитаторах. В этом 
случае постановка задачи выглядела следующим образом: 


д.Т-У.(кУТ)+у.УТ=0, ар 


где Г — температура; { — время; к — коэффициент температуропроводности. 
Начальное распределение температуры: 


Т(Е=0,х)=1 (<), (12) 


было принято следующим: внутри лопатки постоянная температура 400 °С, температура окружающей 
среды 0 °С. Граничные условия отвечали отсутствию теплового потока на внешнем контуре С: 


т 


=0,т 0. (13) 


Супт<0 = 
дп Су-п> 0 


Распределение скоростей потока в задаче (11) вычислялось при решении задачи (10), при этом скорость на 
контуре С принята равной у= 10[Шс] (размерность длины в единицах длины хорды). Коэффициенты 
теплопроводности лопатки и газа приняты равными 0,1 и 0,01 [[/с] соответственно. 

Программный код для решения задач (10), (11) был написан на входном языке универсального конечно- 
элементного пакета ЕгееЕет-+ [20]. 
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На рис.2 приведены рассчитанные в результате решения задач (11), (12) распределения давления и 
температуры при остывании лопатки в течение 25 с. в аэродинамическом потоке для трех изученных профилей 
лопаток, представленных на рис. |1. В таблице2 приведены вычисленные по ним значения минимальных 
относительных давлений вблизи поверхностей лопаток и максимальные температуры лопаток. 


в) 


Рис. 2. Рассчитанные распределения давления (р, слева) и температуры (Т, справа) при моделировании обтекания и 
остывания в потоке лопаток исследованных профилей (обозначения профилей — как на рис. 1). Угол направления вектора 
потока с направлением хорды лопатки составлял 10° 


Таблица 2 
Минимальное относительное давление вблизи поверхности лопатки 
и максимальная температура лопатки после 25 с остывания 
в аэродинамическом потоке 


НО Профиль лопатки 
[13] [18] [19] 
Ри [| Ро 0,599 0,699 0,720 
Тик (#=25), °С 156,9 2452 263,4 


Обсуждение и заключение. Судя по данным, представленных в таблице 1, коэффициенты детерминации В? 
для всех изученных профилей лопаток близки к единице несмотря на то, что формы профилей существенно 
различались. Это подтверждает хорошее качество аппроксимации и иллюстрируется также данными на рис. 1. 
Таким образом, предложенный метод позволяет с высокой точностью построить аппроксимирующие линии 
профилей лопаток для различных типов ГТД. 

В работе также проиллюстрирована возможность использования аппроксимирующих функций при расчете 
распределения давлений и температур на поверхности лопатки на основании решения задачи газовой динамики. 
Далее данное распределение может быть использовано при оптимизации структуры специальных покрытий 
лопатки, например теплозащитных покрытий. Данные, представленные на рис. 2 и в таблице 2, подтверждают, 
что данная задача также может быть успешно решена. Характеры распределений давлений и температур для 
разных профилей лопаток имеют общие черты, но детальная картина различается. Видно в частности, что лопатка 
профиля [19], имеющая большую, по сравнению с остальными, площадь сечения, характеризуется большим 
значением относительного давления в зоне выходной кромки и меньшей скоростью остывания. Это 
иллюстрируется также данными, приведенными в таблице 2. Поскольку с точки зрения прочности адгезионной 
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связи покрытия с подложкой наиболее опасными являются усилия, направленные нормально от поверхности 
лопатки, то актуальным является относительное минимальное давление вблизи поверхности, значения которого 
представлены в таблице 2. На основании полученного расчета можно обоснованно подойти к проектированию 
специального покрытия поверхности лопатки конкретного профиля: регулировать материалы и толщину 
покрытия в зависимости от степени нагруженности отдельных зон поверхности. 

Заметим, что данные результаты расчетов представлены лишь в качестве иллюстрации возможностей 
предложенной методики аппроксимации профилей. Поэтому использованная в работе модель газовой динамики 
была достаточно упрощенной. В конкретной задаче необходимо использовать более строгую модель с 
параметрами, отвечающими заданному типу и режиму работы ГТД. 

Таким образом, предлагаемая методика аппроксимации профиля лопатки ГТД, основанная на использовании 
ортогональных многочленов, является удобным инструментом, позволяющим автоматизировать создание 
геометрической модели лопатки и проводить сравнение различных типов и профилей лопаток, решая 
соответствующие задачи газовой динамики. При этом для заданного профиля лопаток и условий работы ГТД 
можно получить распределение температур и усилий, действующих на поверхности лопатки, что необходимо 
для прогнозирования долговечности специальных покрытий [21]. На основе данного расчета можно 
оптимизировать технологию нанесения покрытия, варьируя его толщину и/или состав, увеличивая коэффициент 
запаса долговечности в наиболее опасных областях. 
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